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Abstract 

We estimate the frequency spectrum of the diffuse anomalous microwave emission 
(AME) on the North Celestial Pole (NCP) region of the sky with the Correlated Com- 
ponent Analysis (CCA) component separation method applied to WMAP 7-yr data. 
The NCP is a suitable region for this analysis because the AME is weakly contami- 
nated by synchrotron and free-free emission. By modeling the AME component as a 
peaked spectrum we estimate the peak frequency to be 21.7 it 0.8 GHz, in agreement 
with previous analyses which favored fp < 23 GHz. The ability of our method to cor- 
rectly recover the position of the peak is verified through simulations. We compare the 
estimated AME spectrum with theoretical spinning dust models to constrain the hy- 
drogen density nu- The best results are obtained with densities around 0.2-0.3 cm~^, 
typical of warm ionised medium (WIM) to warm neutral medium (WNM) conditions. 
The degeneracy with the gas temperature prevents an accurate determination of njj, 
especially for low hydrogen ionization fractions, where densities of a few cm~^ are also 
allowed. 

1 Introduction 



The anomalous microwave emission (AME) component is highly correlated with the far infra- 
red dust emission ([21], |23], [2], |22], [12], [6], [9], [30], [TT]) and is believed to be the 
electric dipole radiation from small spinning dust grains ([lO]). Spinning dust models predict 
a peaked spectrum ranging 10-150 GHz depending on the local physical conditions ([T|.[3T]). 
Probing the peak of the emission enables us to compare models with observations and is 
the best way to distinguish the AME from synchrotron emission and free-free emission, both 
having a power-law spectrum. 

The AME has been studied in individual dust clouds associated with reflection nebulae, 
molecular clouds, photo-dissociation regions and HII regions ([14],[26]). The best examples 
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of peaked AME spectra are probably the Perseus and p Ophiuchi molecular clouds, which 
have been probed with high accuracy thanks to the availability of many different datasets. 
In these regions the derived peak is at ~30 GHz ([26]). 

The study of the AME in more diffuse regions essentially becomes a component separa- 
tion problem, which is often complicated by the lack of low-frequency data covering large 
sky areas. Diffuse dust-correlated emission has been detected with COBE-DMR{[^) and 
WMAP data ([6], [1], [9], [25], [15]), however the details of the spectrum are still unclear. 
Template-fitting analyses show that the dust-correlated emission between 20 and 60 GHz is 
well described by a power-law ([2], [B], [H], [I5])- Given the error bars, a peaked spectrum is 
not ruled out, but a low peak frequency (< 23 GHz) is favored. Alternatively, the observed 
spectrum could result from the superposition of multiple peaked components. The compo- 
nent separation analysis done by [1] found that a peaked AME model gives better results in 
terms of CMB cleaning, however no estimation of the AME peak frequency was attempted. 

In this work we address the estimation of the peak frequency of the diffuse AME on 
WMAP 7-yr data complemented by ancillary data with the Correlated Component Analysis 
(CCA, [27J) component separation method. We consider the North Celestial Pole (NCP) 
region of the sky, centered on Galactic coordinates {l,b) ~ (125°, 25°). As already noted 
by [6], this region of the sky is particularly suited for this analysis since there is significant 
thermal and anomalous dust emission, while synchrotron emission and free-free emission 
(traced by the 408 MHz from [i9\ and the Ha map from [8j) are faint. 

2 Description of the method 

In the following we describe the principles of operation of the harmonic-domain CCA; for 
further details we refer the reader to [27]. The CCA ([3], [27]) exploits second-order statistics 
to estimate the frequency scaling of the components from the statistics of data and noise. 
For each position on the sky (each pixel) we write our data model as: 

x(r) = [B*Hs](r) + n(r), (1) 

where s is a vector whose elements contain the different components (CMB, free-free emis- 
sion, synchrotron emission, thermal dust emission, AME); x and n are vectors too, each 
element containing respectively the data and the instrumental noise for each frequency. The 
diagonal matrix B contains the instrumental beams for all frequency channels and * denotes 
convolution. The matrix H is called the mixing matrix and it contains the frequency scaling 
of the components for all the considered frequencies. 

By translating eq. ([1]) in the harmonic domain, the convolution becomes a multiplication 
and the data model becomes a linear mixture. For each transformed mode (£) we can write: 

X(£) = B(£)HS(£) + N(£), (2) 

where X, S, and N are the transforms of x, s, and n, respectively, and B is the transform 
of the matrix B. The cross-spectra of the data Cx(^), sources Cs(^) and noise, Cn(^), are 
related by: 

Cx(^) - C„(£) = B(£)HCs(£)H^Bt(£) (3) 

where the dagger superscript denotes the adjoint matrix. 

The left-hand side of eq. ([3]) can be estimated from the data for a suitable set of spectral 
bins i] this is used by CCA to estimate the mixing matrix and the source cross-spectra 
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Figure 1: Maps of the region of interest. Top: WMAP K band (left), 94 GHz dust prediction 
(right, |TT]). Bottom: free-free (left) and synchrotron (right) predictions at 23 GHz based on 
the Ha map ([S]) and the 408 MHz map ([H]) assuming a spectral index (3^ = —2.9. 
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on the right-hand side of the equation. Due to the scahng ambiguity of the problem, the 
mixing matrix is normahzed at a reference frequency. To reduce the number of unknowns, 
the mixing matrix is parametrized through a parameter vector p (such that H = H(p)), 
by adopting suitable fitting relations for the spectra of the astrophysical components (as 
detailed in the next section). 

By using 2-dimensional discrete Fourier transforms it is possible to apply the harmonic- 
domain CCA to square sky patches. This approach is advisable as the mixing matrix varies 
on the sky. The HEALPix ([H]) data on the sphere are projected on the plane tangential to 
the center of the patch and re-gridded with a suitable number of bins in order to correctly 
sample the original resolution. The patch size is obtained as a trade-off between the need 
of having uniform spectral properties of the foregrounds and enough statistics for a robust 
computation of the auto- and cross-spectra of the data. In the present analysis we use a 
patch size of 30° x 30°. By using higher resolution data (e.g. Planck data) it is possible to 
reduce the patch size by carrying the analysis up to a larger multipole. We analyzed the 
sky patch centered on [l, b) = (124°, 26°) (shown in Fig. [T]) and verified the stability of the 
results for a sample of patches shifted in latitude and longitude up to 10°. 

3 Description of the analysis 

We used the following datasets: 

• WMAP 7-yr K, Ka, Q, V and W bands ([20]). All maps have been used at the original 
resolution except the K band, which has been smoothed to 1° resolution to reduce the 
effect of beam asymmetry; 

• 408 MHz map ([IH]) to trace the synchrotron component; 

• Predicted free-free emission at 23 GHz based on the Ha map by [8J corrected for dust 
absorption with the E(B-V) map by [2H] assuming = 0.3; 

• Predicted dust emission at 94 GHz by [11] . 

The different resolution of the data maps is accounted for with the beam matrix B in 
eq. (|3]). In particular, the beams are deconvolved for the mixing matrix estimation. The noise 
properties for the WMAP majps have been computed by simulating different noise realizations 
starting from the input Nohs maps and ctq values. For the maps used as foreground templates 
(408 MHz map, 23 GHz free-free map and 94 GHz dust prediction) we assumed a Gaussian 
noise at the 10% level, which is much higher than the instrumental noise. This extra noise 
mimics the error on the template as a tracer of the true component, the 10% level being 
indicative. We verified, however, that the results are not sensitive to the exact value assumed. 
We modeled the data as a mixture of five components: CMB, synchrotron emission, thermal 
dust emission, free-free emission and AME. Each frequency map contributes five elements to 
the mixing matrix, one for each component. The synchrotron, free- free and dust templates 
are used as frequency channels having only one mixing matrix entry, that of the corresponding 
component, different from zero. 

It is worth noting that, compared with previous work, we are not assuming correlation 
between the AME and thermal dust emission and we are not exploiting any template for 
the AME. Though the AME and the thermal dust emission are significantly correlated, the 
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Table 1: Parameters for the spectra in Fig. [2j The SpDust input parameters are: the 
total hydrogen number density nn, the gas temperature T, the intensity of the radiation 
field relative to the average interstellar radiation field x? the hydrogen ionization fraction 
= nii+frin and the ionized carbon fractional abundance xc = nc+Z'^-n- We also report 
the CCA spectral parameters Up, meo [eq. 



model name 


rzH [cm ^ 


T[K] 


X 








meo 


WNM 


0.4 


6000 


1.00 


0.10 


0.0003 


24.34 


7.64 


CNM 


30.0 


100 


1.00 


0.0012 


0.0003 


29.26 


5.31 


WIM 


0.1 


8000 


1.00 


0.99 


0.001 


27.64 


5.99 


MC 


300 


20 


0.01 


0.0 


0.0001 


38.13 


2.23 



correlation is not supposed to be perfect, as, according to spinning dust models, the AME 
traces the distribution of smaller dust grains (PAHs). Indeed, the intensity ratio of AME 
and thermal dust is found to vary in the sky by a factor of ~2 (e.g., [H]). 

We estimated the frequency spectrum of the AME while assuming the spectra of the other 
components to be known. For CMB we used the usual black-body law with temperature 
of 2.726 K ([I3]); for free-free we used the model Tpj,j^fj(i/) cx G x (z//10)~^ where G = 
3.96(T4)°'^^(i//40)~°'^^ is the Gaunt factor, which is responsible for the departure from a 
pure power-law behavior, and T4 = 0.7 is the electron temperature Tg in units of 10^ K. The 
thermal dust emission has been modeled as a grey-body with temperature = 18K and 
spectral index /3d = 1.7. This is consistent with the 1-component dust model by Planck 
data, extending the observations up to 857 GHz, will allow a refinement of this model. For 
synchrotron we assumed a power-law scaling with spectral index Ps = —2.9. 

For the AME we used the parametric relation proposed by [4j: 

logTRj,amc O^ OC + 2 log Z/ + , 4 

Vlog(z/p/60) / 21og(z/p/60) 

which is a parabola in the log(i^) — log(S') plane parametrized in terms of peak frequency Up 
and slope at 60 GHz mgo- The CCA solves for the two parameters meo and z/p, so the accuracy 
of the estimation of the AME peak frequency Up is related to the ability of the parametric 
relation to reproduce the true spectrum. We verified, however, that this parametric model 
can fit a wide range of spinning dust spectra. This is shown in Fig. [21 comparing spinning 
dust spectra produced with the SpDust ([Ij, [29j) code with the results of the fit of eq. (jl]) 
obtained by minimizing the for the set of frequencies considered in this work. The physical 
models that we consider are: warm neutral medium (WNM), cold neutral medium (CNM), 
warm ionised medium (WIM), and molecular cloud (MC). In general, the fits are accurate 
up to 61 GHz, while at 94 GHz the parametric relation may not be able to reproduce the 
input spectra in detail. This is a consequence of fitting complex spectra with only a few 
parameters, the fit being less accurate where the AME is weaker. The best-fit parameters 
that we obtain, reported in Table [H vary significantly from one input model to another. 

4 Test on simulated data 

We simulated WMAP 7-yr data by assuming monochromatic band-passes centered at the 
central frequency of the realistic bands, Gaussian beams at the nominal values, Gaussian 
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Figure 2: Theoretical spinning dust models produced with SpDust (solid lines) and fitted 
with CCA parametric model (triangles, same color of the model). Input SpDust parameters 
and best-fit parameters for the CCA model are provided in Table [H 



noise generated according to realistic (spatially varying) RMS. Our model of the sky consists 
of the following components: 

• CMB emission constrained by the best-fit power spectrum model to WMAP 7-yr; 

• Synchrotron emission given by [19] scaled in frequency with a power-law model with a 
spatially varying synchrotron spectral index /3s as modeled by [IB] ; 

• Free-free component modeled by the Ha map of |8j corrected for dust absorption with 
the E(B-V) map from [28] with /d = 0.3 and scaled in frequency with a power-law 
with fixed spectral index of —2.14; 

• Thermal dust emission modeled with the 100 /xm map from [28] scaled in frequency 
with the best-fit model by [H] which consists of two grey-body laws having different 
temperatures and emissivity indices; 

• AME emission modeled by the E(B— V) map from [28] with intensity at 23 GHz cal- 
ibrated on the results by [6] for average intermediate-latitude conditions and scaled 
according to SpDust models. 

It is worth noting that this sky model is more complex than the spectral model assumed 
in the component separation. The synchrotron spectral index is spatially varying and both 
the thermal dust AME spectra, though spatially constant, are not generated with the CCA 
parametric models. This has been done purposely, to reflect a more realistic situation. To 
test the ability of our pipeline to recover the AME spectrum in different situations, we 
adopted two SpDust models for the AME: one peaking at 26 GHz and the other peaking at 
19 GHz. 
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Figure 3: Comparison of true input AME spectra (black solid lines) with CCA estimates 
(red dashed lines) for two AME models, one peaking at 26 GHz (left) and one at 19 GHz 
(right). The red shaded areas show the recovered z/p and mgo with la error bar. The results 
for the peak frequency are z/p = 26.5 ± 0.5 and Up = 18.5 ± 1.5 for the left and right panels 
respectively. The grey vertical lines indicate the frequency of 23 GHz which sets the limit of 
our data. 



We performed an estimate of the AME spectrum on the simulated data following the same 
procedure applied to the real data. The results are presented in Fig. |3l For AME peaking 
at 26 GHz the recovery of the spectrum is very accurate; the result on the peak frequency 
is Up = 26.5 ± 0.5. The accuracy of the estimation is helped by the weak synchrotron and 
free-free contamination in this region of the sky (in a more general case the error on z/p 
would be of a few GHz). With the AME peaking at 19 GHz the errors increase, but the 
recovery of the spectrum is still satisfactory; we obtain z/p = 18.5 ± 1.5. Our pipehne is able 
to distinguish very clearly between the two input models and to correctly estimate the peak 
frequency of the AME in both cases. 

5 Results 

The results for the AME spectrum on the NCP are showed in Fig. HJ the estimated peak 
frequency is z/p = 21.7±0.8. The diamonds with error bars are the AME mixing matrix 
elements at WMAP frequencies and related uncertainties derived from the results on z/p and 

By comparing these data points with theoretical models produced with SpDust we can 
constrain the physical properties of the medium with the hypothesis of spinning dust emis- 
sion. Some relevant parameters of the SpDust code are: hydrogen number density nn (cm~^); 
gas temperature T (K); intensity of the radiation field with respect to average Galactic con- 
ditions, X; and hydrogen ionization fraction, xh- Advanced parameters describe the dust 
grains and define the abundance of other elements. 
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Figure 4: Best-fit AME model (black line) and results in terms of Up and mgo considering 
Icr uncertainties (red areas). The peak frequency results Up = 21.7 ± 0.8. Diamonds with 
error bars are the results in terms of normalized intensity at WMAP bands and related 
uncertainties. 



The hydrogen density nn is sensitive to the position of the peak and in general to the 
shape of the spectrum in the considered frequency range. It is also degenerate with other pa- 
rameters, in particular the gas temperature T. We investigated this degeneracy by sampling 
the 2-dimensional rin-T space with a grid approach while conditioning the remaining param- 
eters. For each point in the parameter space we computed the likelihood as A = exp(— x^/2) 
where is the standard chi-square. In computing the we did not consider the 94 GHz 
data point, which could be biased, as verified through simulations (see Sect. H]). We repeated 
the likelihood estimation for different choices of the parameters xh and x (varying all those 
parameters simultaneously would be computationally too demanding). As the CCA outputs 
are in normalized units, we compared both the model and the data after scaling each of them 
by the sum of the intensities at 23, 33, 41 and 61 GHz. 

The results of the nn-T parameter estimation are shown in Fig. [5l The colored lines show 
the likelihoods for tt-h for different values of T as detailed in the legend, while the black lines 
show the likelihood for rin marginalized over T in the range 100-10000 K. Different panels 
correspond to different choices for the other parameters: xh = 0.001, x = on the top left; 
xh = 0.001, X = 0.1, on the top right, xh = 0.01, x = on the bottom left and xh = 0.1, 
X = 1 on the bottom right. The marginalized statistics for n-n are reported in Table [2] (An^ 
and Any being the lower and upper la error on nn). The uyi-T degeneracy weakens the 
constraint on rin in particular on the high-density side of the distribution. The degeneracy 
decreases with increasing hydrogen ionization fraction. 

The best results are obtained for densities of 0.2-0.4 cm~^, typical of WIM/WNM condi- 
tions. In Fig. [6] we show the data compared to some SpDust models to exemplify the quality 
of the fit. The comparison of the red (tt-h = 0.25 cm~'^, T = 6000 K, xh = 0.1, x = 1) and 
magenta (nn = 0.25 cm~^, T = 6000 K, xh = 0.01, x = 1) lines illustrates the effect of xh- 
The blue line corresponds to riu = 0.3 cm~^, T = 6000 K, xh = 0.31, x = 1; compared to the 
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Figure 5: Colored lines: likelihoods for nn for different values of T as detailed in the legend; 
black line: likelihood for nn marginalized over T. Top left: Xh = 0.001, x = 1; top right: 
same of top left with x = 0.1; bottom: same as top left with xh = 0.01 (left) and xn = 0.1 
(right). 
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Figure 6: Comparison between data (diamonds and error bars) and different SpDust models. 
Red: = 0.25 cm~^ T = 6000 K, = 0.1, x = 1; magenta: same as red with xu = 0.01; 
blue: riH = 0.30 cm'^ T = 6000 K, xn = 0.01, x = 1; green: nu = 0.40 cm-^ T = 6000 K, 
xh = 0.001, X = 1; cyan: same as blue with x = O-l- 



Table 2: Statistics for nn marginalized over T corresponding to the results of FigJS] 



xn 


X 


nn 






0.001 


1 


0.3 


0.1 


0.9 


0.01 


1 


0.3 


0.1 


0.7 


0.1 


1 


0.25 


0.09 


0.07 


0.001 


0.1 


0.4 


0.1 


1.2 



magenta line, it illustrates the effect of increasing riH- Finally, the cyan (nn = 0.40 cm ^, 
T = 6000 K, xh = 0.001, x = 0.1) and green (nn = 0.40 cm'^, T = 6000 K, xh = 0.001, 
X = 1) lines show the effect of x- The availability of data at frequencies of a few GHz will 
be very valuable to improve the constraints on such parameters. 

To complement the physical description of the medium we computed the hydrogen column 
density (cm~^) from the brightness at 100 yUm, Jioo, through Jioo [MJysr~^] = (0.69 ± 
0.03) X A^H [10^° cm-2] ([S]). We used the IRIS ([24]) band 4 (100 /im) data and integrated the 
emission inside a 1° beam centered in I = 25°, b = 125°, corresponding to bright dust emission 
(see Fig. [1]), to compute a representative value for Jioo- This analysis yields A^h = (2.9 ± 
0.2) X 10^°cm~^, which is a plausible value for the diffuse ISM environment at intermediate 
latitudes. 



6 Conclusions 



We have studied the spectrum of the diffuse AME with WMAP 7-yr and ancillary data in 
the North Celestial Pole (NCP) region of the sky. In this region the AME dominates the low- 
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frequency emission, as both synchrotron and free- free are faint. Previous template-fitting 
analysis by [6J found that the AME spectrum in this region is consistent with a power-law; 
the same would apply to most diffuse AME at intermediate latitudes ([2], [6], [4], [9], [25], 
[15j). This favors a low peak frequency (z/p < 23 GHz). 

For our analysis we rely on the CCA component separation method, which exploits the 
data auto- and cross-spectra to estimate the frequency spectra of the components in terms 
of a set of spectral parameters. Our method models the AME as a peaked spectrum and fits 
for the peak frequency, Up, and the slope at 60 GHz, mgo. We verified with simulations that 
we are able to correctly recover the AME spectrum and more specifically the peak frequency, 
even when it is below 23 GHz. We get z/p = 21.7 ± 0.8 GHz, which is both a confirmation 
and an improvement of the previous results. This result relies on the assumption that our 
parametric model for the AME is a good representation of the true spectrum, which has 
been verified for a wide range of theoretical spinning dust models. 

Using the SpDust code, we linked the estimated spectrum to the local physical conditions 
with the hypothesis of spinning dust emission. We investigated in particular the hydrogen 
density (nn [cm~'^]), which is sensitive to the position of the peak, and its degeneracy with 
the gas temperature T [K]. The densities that we get are those typical of WIM/WNM con- 
ditions (0.2-0.4 cm~^). For a low hydrogen ionization fraction (nn < 10~^), densities up to 
a few cm~^ are allowed by the nn-T degeneracy. Lower radiation fields (x) require slightly 
lower nn. By considering 100 /xm data for our region we obtain a hydrogen column density 
A'h ~ lO^'^cm"^. Overall, the recovered AME spectrum is found to be consistent with that 
predicted by spinning dust models for plausible physical conditions. 

7 Acknowledgements 

S. Ricciardi acknowledges support by MIUR through PRIN 2009 grant n. 2009XZ54H2 

References 

[1] Y. Ali-Hai'moud, C. M. Hirata, and C. Dickinson. A refined model for spinning dust 
radiation. MNRAS, 395:1055-1078, May 2009. 

[2] A. J. Banday, C. Dickinson, R. D. Davies, R. J. Davis, and K. M. Gorski. Reappraising 
foreground contamination in the COBE-DMR data. MNRAS, 345:897-911, Nov. 2003. 

[3] A. Bonaldi, L. Bedini, E. Salerno, C. Baccigalupi, and G. de Zotti. Estimating the 
spectral indices of correlated astrophysical foregrounds by a second-order statistical 
approach. MNRAS, 373:271-279, Nov. 2006. 

[4] A. Bonaldi, S. Ricciardi, S. Leach, F. Stivoli, C. Baccigalupi, and G. de Zotti. WMAP 
3-yr data with Correlated Component Analysis: anomalous emission and impact of 
component separation on the CMB power spectrum. MNRAS, 382:1791-1803, Dec. 
2007. 

[5] F. Boulanger, A. Abergel, J. -P. Bernard, W. B. Burton, F.-X. Desert, D. Hartmann, 
G. Lagache, and J.-L. Puget. The dust/gas correlation at high Galactic latitude. A&A, 
312:256-262, Aug. 1996. 



Spectrum of the Anomalous Microwave Emission in the North Celestial Pole with WMAP 7-Year datal2 



[6] R. D. Davies, C. Dickinson, A. J. Banday, T. R. Jaffe, K. M. Gorski, and R. J. Davis. 

A determination of the spectra of Galactic components observed by the Wilkinson Mi- 
crowave Anisotropy Probe. MNRAS, 370:1125-1139, Aug. 2006. 

[7] A. de Oliveira-Gosta, M. Tegmark, R. D. Davies, G. M. Gutierrez, A. N. Lasenby, 
R. Rebolo, and R. A. Watson. The Quest for Microwave Foreground X. ApJ, 606:L89- 
L92, May 2004. 

[8] G. Dickinson, R. D. Davies, and R. J. Davis. Towards a free-free template for GMB 
foregrounds. MNRAS, 341:369-384, May 2003. 

[9] G. Dobler and D. P. Finkbeiner. Extended Anomalous Foreground Emission in the 
WMAP Three- Year Data. ApJ, 680:1222-1234, June 2008. 

[10] B. T. Draine and A. Lazarian. Electric Dipole Radiation from Spinning Dust Grains. 
ApJ, 508:157-179, Nov. 1998. 

[11] D. P. Finkbeiner, M. Davis, and D. J. Schlegel. Extrapolation of Galactic Dust Emission 
at 100 Microns to Gosmic Microwave Background Radiation Frequencies Using FIRAS. 
ApJ, 524:867-886, Oct. 1999. 

[12] D. P. Finkbeiner, G. I. Langston, and A. H. Minter. Microwave Interstellar Medium 
Emission in the Green Bank Galactic Plane Survey: Evidence for Spinning Dust. ApJ, 
617:350-359, Dec. 2004. 

[13] D. J. Fixsen, J. L. Weiland, S. Brodd, M. G. Hauser, T. Kelsall, D. T. Leisawitz, J. C. 
Mather, K. A. Jensen, R. A. Schafer, and R. F. Silverberg. Gomparison of the GOBE 
FIRAS and DIRBE Galibrations. ApJ, 490:482-+, Dec. 1997. 

[14] R. Genova-Santos, R. Rebolo, J. A. Rubino- Martin, G. H. Lopez-Garaballo, and S. R. 
Hildebrandt. Detection of Anomalous Microwave Emission in the Pleiades Reflection 
Nebula with Wilkinson Microwave Anisotropy Probe and the GOSMOSOMAS Experi- 
ment. ApJ, 743:67, Dec. 2011. 

[15] T. Ghosh, A. J. Banday, T. Jaffe, G. Dickinson, R. Davies, R. Davis, and K. Gorski. 
Foreground analysis using cross-correlations of external templates on the 7-ycar Wilkin- 
son Microwave Anisotropy Probe data. MNRAS, 422:3617-3642, June 2012. 

[16] G. Giardino, A. J. Banday, K. M. Gorski, K. Bennett, J. L. Jonas, and J. Tauber. 
Towards a model of full-sky Galactic synchrotron intensity and linear polarisation: A 
re-analysis of the Parkes data. A&A, 387:82-97, May 2002. 

[17] B. Gold, N. Odegard, J. L. Weiland, R. S. Hill, A. Kogut, G. L. Bennett, G. Hinshaw, 
X. Chen, J. Dunkley, M. Halpern, N. Jarosik, E. Komatsu, D. Larson, M. Limon, S. S. 
Meyer, M. R. Nolta, L. Page, K. M. Smith, D. N. Spergel, G. S. Tucker, E. Wollack, 
and E. L. Wright. Seven-year Wilkinson Microwave Anisotropy Probe (WMAP) Obser- 
vations: Galactic Foreground Emission. ApJS, 192:15, Feb. 2011. 

[18] K. M. Gorski, E. Hivon, A. J. Banday, B. D. Wandelt, F. K. Hansen, M. Reinecke, and 
M. Bartelmann. HEALPix: A Framework for High-Resolution Discretization and Fast 
Analysis of Data Distributed on the Sphere. ApJ, 622:759-771, Apr. 2005. 



Spectrum of the Anomalous Microwave Emission in the North Celestial Pole with WMAP 7-Year datalS 



[19] C. G. T. Haslam, C. J. Salter, H. Stoffel, and W. E. Wilson. A 408 MHz all-sky 
continuum survey. II - The atlas of contour maps. A&AS, 47:1, Jan. 1982. 

[20] N. Jarosik, C. L. Bennett, J. Dunkley, B. Gold, M. R. Greason, M. Halpern, R. S. 
Hill, G. Hinshaw, A. Kogut, E. Komatsu, D. Larson, M. Limon, S. S. Meyer, M. R. 
Nolta, N. Odegard, L. Page, K. M. Smith, D. N. Spergel, G. S. Tucker, J. L. Weiland, 
E. Wollack, and E. L. Wright. Seven-year Wilkinson Microwave Anisotropy Probe 
(WMAP) Observations: Sky Maps, Systematic Errors, and Basic Results. ApJS, 192:14, 
Feb. 2011. 

[21] A. Kogut. Spatial Correlation Between Dust and Diffuse Ionized Gas. In American 
Astronomical Society Meeting Abstracts, volume 28 of Bulletin of the American Astro- 
nomical Society, page 1295, Dec. 1996. 

[22] G. Lagache. The large-scale anomalous microwave emission revisited by WMAP. A&A, 
405:813-819, July 2003. 

[23] E. M. Leitch, A. C. S. Readhead, T. J. Pearson, and S. T. Myers. An Anomalous 
Component of Galactic Emission. ApJ, 486:L23, Sept. 1997. 

[24] M.-A. Miville-Deschenes and G. Lagache. IRIS: A New Generation of IRAS Maps. 
In L. Armus & W. T. Reach, editor. Astronomical Society of the Pacific Conference 
Series, volume 357 of Astronomical Society of the Pacific Conference Series, page 167, 
Dec. 2006. 

[25] M.-A. Miville-Deschenes, N. Ysard, A. Lavabre, N. Ponthieu, J. F. Macias-Perez, J. Au- 
mont, and J. P. Bernard. Separation of anomalous and synchrotron emissions using 
WMAP polarization data. A&A, 490:1093-1102, Nov. 2008. 

[26] Planck Collaboration. Planck early results 20: New light on anomalous microwave 
emission from spinning dust grains. Submitted to A&A, [arXiv:astro-ph/1101.2031], 
2011. 

[27] S. Ricciardi, A. Bonaldi, P. Natoli, G. Polenta, C. Baccigalupi, E. Salerno, K. Kayabol, 
L. Bedini, and G. de Zotti. Correlated component analysis for diffuse component sepa- 
ration with error estimation on simulated Planck polarization data. MNRAS, 406:1644- 
1658, Aug. 2010. 

[28] D. J. Schlegel, D. P. Finkbeiner, and M. Davis. Maps of Dust Infrared Emission for Use 
in Estimation of Reddening and Cosmic Microwave Background Radiation Foregrounds. 
ApJ, 500:525-+, June 1998. 

[29] K. Silsbee, Y. Ali-Hai'moud, and C. M. Hirata. Spinning dust emission: the effect of 
rotation around a non-principal axis. MNRAS, 411:2750-2769, Mar. 2011. 

[30] N. Ysard, M. A. Miville-Deschenes, and L. Verstraete. Probing the origin of the mi- 
crowave anomalous foreground. A&A, 509:L1, Jan. 2010. 

[31] N. Ysard and L. Verstraete. The long-wavelength emission of interstellar PAHs: char- 
acterizing the spinning dust contribution. A&A, 509:A12, Jan. 2010. 



